Orbital Order and Spontaneous Orthorhombicity in Iron Pnictides 
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A growing list of experiments show orthorhombic electronic anisotropy in the iron pnictides, in 
some cases at temperatures well above the spin density wave transition. These experiments include 
neutron scattering, resistivity and magnetoresistance measurements, and a variety of spectroscopies. 
We explore the idea that these anisotropies stem from a common underlying cause: orbital order 
manifest in an unequal occupation of dxz and dyz orbitals, arising from the coupled spin-orbital 
degrees of freedom. We emphasize the distinction between the total orbital occupation (the inte- 
grated density of states), where the order parameter may be small, and the orbital polarization 
near the Fermi level which can be more pronounced. We also discuss light-polarization studies of 
angle-resolved photoemission, and demonstrate how x-ray absorption linear dichroism may be used 
as a method to detect an orbital order parameter. 

PACS numbers: 72.10.Fk, 74.70.Dd, 75.25.Dk, 78.70.Dm 



The iron pnictides"'^ possess a phase diagram similar to 
the cuprates with superconductivity in close proximity 
to magnetic order. Upon doping the magnetism is sup- 
pressed and superconductivity emerges.^ In contrast, the 
pnictide parent compounds are metals, rather than insu- 
lators, with a Fermi surface consisting of bands mostly 
from the Fe id t2g orbitals. These orbitals form two hole 
pockets with mostly d^z and dyz characters around the 
Brillouin zone (BZ) center, and two electron pockets with 
additional d^y weight at the zone corners.'^ 

Whether a weak-coupling itinerant electron or a 
strong-coupling localized moment approach better de- 
scribes the pnictides remains a topic of debate. It is gen- 
erally accepted that the Fe on-site, intra-orbital interac- 
tion f/fc^'^ is smaller than the overall Fe 3d bandwidth ^ 
4eV. In particular, x-ray absorption spectroscopy (XAS) 
indicates U ~ 2eV.^ In this regard, the pnictides are 
weakly, or at most moderately, correlated materials. On 
the other hand, both itinerant and local pictures seem 
to well describe the neutron scattering measurements.^''' 
The above observations therefore suggest a coexistence 
of localized moments and itinerant electrons in the iron- 
pnictide materials.^ They have aspects such as metal- 
licity where correlations play a minor role, while anti- 
ferromagnetism and local properties derive directly from 
a strong Hund's coupling.^ This dual character is also 
implied by other experiments: Raman scattering finds 
that the dxz/dyz orbitals are more incoherent than dxyf 
de Haas-van Alphen measurements of superconducting 
LaFePO^*^ and the closely related non-superconducting 
compound SrFe2P2^^ reveal that the electron pockets 
have a larger mean free path, possibly tied to the dxy 
band; higher electron mobility is further supported by 
the Hall effect, as well as nuclear magnetic relaxation. 

In addition to this dichotomy, a series of experiments 
highlight the importance of orthorhombic anisotropy in 



these materials. ^'^^"^^ In particular, neutron scatter- 
ing indicates extremely anisotropic couplings between 
the Fe moments;^ scanning tunneling microscopy reveals 
unidirectional electronic nano-structures,^^ and nuclear 
quadrupole resonance shows a local electronic order in 
the Fe layers. These observations of broken C4 tetrago- 
nal symmetry could stem from an unequal occupation of 
the dxz/dyz orbitals, arising from the coupled spin-orbital 
physics. -'^^"^■^ A nematic state is formed once the dxz/dyz 
orbitals are ordered in a translationally invariant way.^^ 
On the other hand, pure spin physics can lead to ne- 
maticity^^ resulting from order by disorder. The latter 
scenario can cause an orbital polarization upon breaking 
C4 symmetry. On symmetry grounds the two scenarios 
are not distinct since both lead to an Ising order parame- 
ter. However, they can be distinguished by the spin-spin 
correlations and spin dynamics above the spin density 
wave (SDW) transition. 

In this study, we assume that orbital occupation is 
a key driving variable, and investigate the experimen- 
tal consequences of a net orbital polarization. In par- 
ticular, we show that recent resistivity anisotropy mea- 
surements^^'^'' could result from an unequal dxz/dyz or- 
bital population. Moreover, the Fermi surface recon- 
struction in angle-resolved photoemission (ARPES)^^''^" 
can be better understood with a (partial) orbital order- 
ing. We also calculate XAS linear dichroism which could 
be used to probe the orbital order parameter. 

In the strong coupling limit, orbital order is easier to 
understand if the dxz and dyz bands are j filled, either 
with electrons or holes. In this case, the local occupation 
of the orbitals can be modelled by an Ising variable.^" 
Such a variable carries a full i?ln2 entropy, and the 
orbital order can be robust and easily detectable. In this 
case, metallic behavior must be due to other orbitals such 
as dxy The order can survive at intermediate couplings. 
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but is reduced due to d-electron itinerancy. 

A more subtle orbital order can arise at intermediate 
couplings in a system where the occupancy of dxz and 
dyz orbitals is closer to half filling. In that case, the two 
orbitals have occupancies of 1 + e and 1 — e. For small e, 
this can no longer be an atomic scale fluctuating variable 
but rather one involving a cluster of at least 1/e atoms. 
Orbital order can arise from an instability that promotes 
an electron in a local region from the dxz to the dyz or- 
bital (or vice versa) . The deviation from unit occupancy 
is distributed throughout the system. Near half filling, 
such an electron transfer would raise the system energy 
due to intra-atomic (Hund's) exchange. A competition 
therefore results between Hund's energy which favors sin- 
gle occupancy of both orbitals with no orbital order, and 
inter-atomic exchange which can support orbital order. 

If the system is initially frustrated, e.g. a J1-J2 model 
with J2 comparable to Ji , this instability to inter-orbital 
electron transfer can be the system's way of reducing 
frustration and lowering the inter-atomic exchange en- 
ergy. This phenomenon also happens in triangular lat- 
tice systems which can orbitally order to avoid frustra- 
tion. '^^ In the pnictides, electron transfers can locally 
break the tetragonal symmetry and couple to lattice dis- 
tortions. This could make the nearest-neighbor exchange 
anti-ferromagnetic (AF) in one direction and ferromag- 
netic (FM) in the other. Thus, collinear AF order will 
simultaneously minimize all the inter-atomic exchanges. 
Lowering of this energy would involve every atom in the 
cluster. Therefore, even if Hund's energy is an order of 
magnitude larger than the inter-atomic exchange energy, 
the latter can dominate in substantially large clusters. 
Apart from anisotropic exchange couplings, an orbital 
polarization may in turn cause other orthorhombic phe- 
nomena, which we discuss below. 

In Ba(Fei_2;Co3;)2As2, a striking in-plane (ab) resis- 
tivity anisotropy Ap/p = {pb - Pa)/[^{pb + Pa)], with 
Pb > Pa, onsets at temperatures well above the struc- 
tural Ts and magnetic Tn ordering temperatures. ^^'^^ 
Below Tjv, the magnetic order is AF along the longer 
a-axis and FM along the shorter &-axis. One thus ex- 
pects a higher resistivity along a-axis due to reduced 
hopping. Moreover, if one considers scattering on spin 
fluctuations, ^^'^^''^^ stronger scattering along the AF a- 
axis with large Q momentum transfer is expected; the 
FM order only gives rise to forward scattering. Both 
considerations should lead to pb < Pa, in contradiction 
to experiment. 

It is also puzzling that Ap increases with decreasing 
T. For scattering on magnetic or lattice fluctuations, Ap 
should weaken at lower T because these bosonic degrees 
of freedom freeze-out. The scattering rate F due to spin 
fluctuations is related to the imaginary part of the spin 
susceptibility, F ^ 9^ J2ciX" where g is the cou- 

pling between itinerant electrons and local spins. The 
divergence of x" at the AF wave vector q = Q may ac- 
count for the increased scattering at T ^ TV,^^ but not 
for T well above or below T^. 
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FIG. 1: (Color Online) (a) Resistivity anisotropy Ap/p and 
(b) orbital polarization {5n) as a function of the effective en- 
ergy splitting Acff between the d^c^/dy^ orbitals calculated 
from the 5-band model, (c) Anisotropy in the density of states 
(DOS) with Aoff = 0.1 eV; the arrow indicates the Fermi level 
at zero doping a; = 0. (d) Ap/p and {5n) versus the chemical 
potential p. Depending on the value of p, a substantial Ap 
could result from a modest {5n). 



We use a 5-band model^"* and investigate whether 
simple impurity scattering, combined with weak orbital 
polarization, can account for the observed resistivity 
anisotropy onset at T > Ts,Tm. In particular, we show 
that both the sign and T dependence of Ap can be qual- 
itatively explained in that way. We study the effect of 
a small orbital polarization {5n) = {uxz — nyz) /[\{nxz + 
nyz)], where Uxz {nyz) is the dxz {dyz) orbital occupa- 
tion. {5n) ^ is mimicked by an effective on-site energy 
splitting: txz = eo + Aoff/2, and eyz = eo - Aoff/2. A 
nonzero {5n) could result from a spontaneous symme- 
try breaking mechanism that onsets at T > Ts,Tn,^^ 
or simply from a small external field X such as applied 
stress^^ in the presence of a finite orbital susceptibility 
Xorb oc d{5n)/dX. As T — > T/v, {5n) and its effect can 
be enhanced by the coupled spin-orbital physics. 

We consider scattering on J-function nonmagnetic im- 
purities of strength W , and calculate the resistivity along 
the a, &-axes using the Kubo formula in the first Born 
approximation,^^ where Ap is weakly dependent on W . 
The potential allows electrons to scatter between differ- 
ent bands but only when the initial and final states have 
the same orbital character. Although thermal fluctu- 
ations arc not introduced explicitly, their main effect is 
expected to reduce {5n) with increasing T, i.e. Acff(r) 
should increase as T decreases. 

In Fig. 1(a), a positive splitting Agff > gives rise 
to the experimentally observed pb > Pa- Figure 1(b) 
displays the corresponding value of {5n) . Aoff > corre- 
sponds to {5n) < (or Uxz < n-yz). On the other hand, 
at the chemical potential {p) the dxz density of states is 
larger than the dyz one: Dxz{p) > Dyz{p) [Fig. 1(c)], 
which is the reason for Ap > 0. In the calculation shift- 



3 



ing II to lower energies would cause sizeable changes in 
SD{fi) = D^^dj,) — Dyz{ii) due to the sharp features in 
the density of states [Fig. 1(c)], while only changing the 
overall {5n) by a small amount. This suggests that a sub- 
stantial Ap could result from a modest orbital polariza- 
tion depending on /i [Fig. 1(d)]. Fermi-surface sensitive 
probes such as transport and ARPES depend strongly on 
the value of SD{e) near the Fermi level (Ep), which can 
be appreciable without a large (5n). 

Mean-field calculations of down-folded Hamiltonians'^^ 
find that the coUincar AF order is accompanied by a non- 
zero (Sn) . A robust orbital order usually occurs when the 
dxz (dyz) occupation is away from half filling. Moreover, 
{5n) depends strongly on interactions. Static Coulomb 
interaction studies suggest that the electron wavefunc- 
tions are ordered along the longer, AF a-axis.^^'^^ On the 
other hand, LDA calculations based on Wannier orbitals 
predict an opposite trend. ^'^ Recent double exchange cal- 
culations find that the wavefunction lobes point in the 
short bond, FM direction due to kinetic energy stabi- 
lization.^'* These results indicate that the sign and the 
magnitude of {Sn) could depend subtly on the interplay 
between the band-structure and interactions. 

In our approach, a nonzero Ap originates from an 
orbital-order induced effective mass anisotropy, which if 
it exists could be detectable by optical conductivity mea- 
surements. We note that at low temperatures the resis- 
tivity anisotropy can be enhanced due to an anisotropy 
in carrier density. Recent a6-plane optical measurements 
have reported a depletion of spectral weight upon en- 
tering the magnetically ordered state for light polariza- 
tion only along the 6-axis, with an enhancement of the 
metallic nature of the charge dynamics related to the 
a-axis response. These observations are indeed consis- 
tent with the dc transport measurements. If one also 
considers scattering on spin fiuctuations at T ^ Ta?,^^ 
orbital order can further affect Ap. If dxy electrons are 
scattered by local moments from the dxz/dyz orbitals, an 
anisotropy in the spin-fermion coupling gyy > Qxx could 
occur, resulting in pi, > Pa- 

Recent low T ARPES measurements suggest a domi- 
nant dxz weight on the hole pockets at the F point. 
Experimentally, the orbital characters can be selected by 
dipole matrix elements with different light polarizations. 
Below we use the 5-band model to qualitatively illustrate 
that a weak SDW field consistent with experiments^" will 
not result in a significant difference in the dxz and dyz 
orbital spectral weight. However, a substantial orbital 
weight reconstruction can be obtained when the SDW 
field is accompanied by orbital ordering. The calcula- 
tions are performed at the mean field level. 

Figure 2(a) shows the spectral function A(k, w = Ep) 
in the paramagnetic phase; the real space unit cell con- 
tains two Fe due to the staggered pnictogen height. Fig- 
ure 2(b) shows A(k, uj = Ep) calculated with a symmetry 
breaking SDW field (=50 meV)."'" The real space unit cell 
now contains 4 Fe, but the shadow bands give rise to only 
weak spectral intensity. Similarly, A(k,uj = Ep) in the 
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FIG. 2: (Color Online) (a) A(k,a; = Ep) in the paramagnetic 
(PM) phase plotted in the one Fe unit cell BZ; the green box 
indicates the true two Fe unit cell BZ. (b) A(k,aj = Ep) in 
the SDW phase; the blue box indicates the true four Fe unit 
cell BZ. The shadow bands give rise to only weak spectral 
weight, (c) A(k,iu — Ep) in the SDW phase with orbital 
ordering (Aefi — 0.1 eV). (d)-(f) The d:^z projected spectral 
functions with SDW and various strength of orbital ordering. 
(g)-(i) Similar plots for the dy^ orbital. When the SDW field 
is accompanied by a partial orbital ordering, the outer (inner) 
pocket is dominated by dxz (dyz) orbitals. All the plots are 
broadened with a 10 meV Lorentzian. 



presence of SDW and a partial orbital order introduced 
by an effective splitting Acff = 0.1 eV is shown in Fig. 
2(c). Figure 2 (d)-(f) arc the dxz projected spectral func- 
tions with the SDW field and various strengths of orbital 
ordering. Figures 2 (g)-(i) are similar plots for dyz- 

A weak SDW field alone does not result in a signifi- 
cant difference in the spectral weight of the dxz and dyz 
orbitals. We note that a simple SDW folding cannot 
account for the ARPES experiments cither.**^ However, 
when the SDW field is accompanied by orbital ordering, 
the outer (inner) pocket is dominated by dxz (dyz) char- 
acter. These results imply an orbital weight reconstruc- 
tion across the Fermi surface, emphasizing the potential 
role of the orbital degrees of freedom. 

Conclusive evidence for orbital ordering in the iron 
pnictides could be obtained by performing x-ray dichro- 
ism measurements. When the system is structurally dis- 
torted or orbitally ordered, the absorption spectrum for 
incident light polarized along the x-axis of the crystal 
frame would be different from that along the y-axis. The 
resultant difference in the two XAS spectra Ix{tj) — lyiu) 
is referred to as the linear dichroism (LD). A temperature 
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FIG. 3: (Color Online) (a) The Fe Ls-edge XAS for x and y 
polarization from FEFF calculation for BaFe2As2 in the low 
temperature T = 20-?^ orthorhombic phase. The LD signal 
(7a; — ly) is multiplied by a factor of ten. (b) Atomic multiplet 
calculation with solely orbital ordering Aofi = 0.1 eV for La- 
edge XAS. In both cases, a ~ 3% LD signal results. The 
structure distortion would saturate at low T, but the orbital 
polarization is expected to increase with decreasing T. 



evolution of the LD signals would be highly suggestive of 
the orbital-ordering hypothesis. 

We focus on Fe L-edge XAS where the Fe 2p core elec- 
trons are excited to unoccupied 3d levels. A recent study 
featured the FEFF software*^ to simulate the polarization 
averaged Fe L-edge XAS in a variety of iron-pnictidc 
materials. Results have been obtained in good agree- 
ment with experiments.^ Similar calculations including 
the XAS polarization dependence are shown in Fig. 3(a). 
The LD signal is ~3% based solely on the orthorhombic- 
ity of the low temperature T = 2QK structure. ^'^ 

As the FEFF code has no direct control over orbital oc- 
cupation, we perform an atomic multiplet calculation to 
simulate the effect of solely orbital ordering. Figure 3(b) 
shows the XAS and a ^^3% LD signal obtained from this 
method with Aoff = 0.1 eV. The LD magnitude would 
increase with a stronger Aoff. We note that the sign and 
the magnitude of the LD signal and {5n) are closely re- 
lated, which may either cancel or enhance the LD signal 



from the structural distortion. 

The LD measurements are complicated by the re- 
quirement of single crystals with de-twinned domains, 
which could be achieved by magnetic^^ or mechanical^^'^^ 
means. Moreover, below Tjv the AF phase also con- 
tributes to the LD spectra, so it is important to distin- 
guish these different contributions. A careful tempera- 
ture dependence should be taken near and across differ- 
ent transition lines. An evolution of the dichroic signal 
could be observed at high temperatures T > Ts,Tm if 
orbital ordering exists. At low T the structure distortion 
saturates,*^ but the orbital polarization is expected to 
increase with decreasing T.^^ Therefore, these different 
effects can potentially be distinguished by the LD mag- 
nitude and line shape. 

In conclusion, we have addressed the in-plane or- 
thorhombic anisotropics in the iron pnictides arising from 
orbital order. We demonstrated that an unequal dxz/dyz 
orbital population can lead to a resistivity anisotropy 
similar to experiments. We also showed that the Fermi 
surface reconstruction in ARPES is better understood 
when the SDW field is accompanied by an orbital polar- 
ization. We made predictions for XAS linear dichroism 
which potentially can probe the orbital order parameter. 
In our view, orbital ordering in these materials arises 
from coupled spin-orbital fluctuations, where the orbital 
variable plays a crucial role. Experiments above the SDW 
transition and theories that address the interplay of spin 
and orbital fluctuations'*^ can further clarify the relative 
importance of the orbital degrees of freedom. 
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